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We analyze complete spectra of the lattice Dirac operator in SU(2) gauge theory and demonstrate that the 
distribution of low-lying eigenvalues is described by random matrix theory. We present possible practical ap- 
plications of this random-matrix universality. In particular, reliable extrapolations of lattice gauge data to the 
thermodynamic limit are discussed. 



1. Introduction 

Based on an analysis of sum rules derived by 
Leutwyler and Smilga |Q for inverse powers of 
the eigenvalues of the QCD Dirac operator in a 
finite volume, Shuryak and Verbaarschot con- 
jectured that the so-called microscopic spectral 
density of the Dirac operator, p s {z), is a univer- 
sal function which can be computed in random 
matrix theory (RMT). This quantity is defined 
by 



(1) 



where p(A) = (J2 n ^(^ ~ ^u))a is the spectral 
density of the Dirac operator averaged over all 
gauge field configurations A, V is the space-time 
volume, and £ is the absolute value of the chiral 
condensate. The chiral condensate can be deter- 
mined using the Banks-Casher relation j| 
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The definition of Eq. ([jj) amounts to a magnifi- 
cation of the region of low-lying eigenvalues by a 
factor of the volume. Since the spacing of small 
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eigenvalues is ~ 1/(UE), this definition leads to 
the resolution of individual eigenvalues. 

A number of studies have presented evidence 
supporting the universality conjecture for p s . 
They are summarized in Ref. . Here, we present 
and analyze the results of a high-statistics lat- 
tice study in SU(2) gauge theory with staggered 
fcrmions which confirms this conjecture in a par- 
ticularly direct way. Most of the details omitted 
here, in particular regarding the lattice calcula- 
tions, can be found in Ref. Q, in the talk by 
S. Meyer, and in the poster by M.E. Berbenni- 
Bitsch, respectively (these proceedings). 

2. Microscopic universality 

Depending on the number of colors and on the 
representation of the fcrmions, one must distin- 
guish three different universality classes which 
were classified in Ref. || . These correspond to the 
three chiral ensembles of RMT, the chiral Gaus- 
sian orthogonal (chGOE), unitary (chGUE), and 
symplectic (chGSE) ensemble. Analytic random- 
matrix results have been obtained for all three 
ensembles for a variety of quantities including 
the microscopic spectral density, the distribu- 
tion of the smallest eigenvalue, and microscopic 
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Figure 1. The distribution of the smallest eigenvalue P(A m in) and microscopic spectral density p s (z) (for 
a 10 4 lattice) and the microscopic spectral two-point function p s (x,y) (for an 8 4 lattice.) The histograms 
represent lattice data; the dashed lines are analytical predictions from random matrix theory. 



spectral correlators. We have worked with stag- 
gered fcrmions in SU(2) for which the appropriate 
RMT ensemble is the chGSE. RMT results for the 
chGSE can be obtained by a slight modification 
of the analytic results obtained by Forrester Q 
and by Nagao and Forrester [jTj for generalized 
Laguerre ensembles. In particular, we have 
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for the microscopic spectral density and 

^(Amin) - y|c( C A min ) 3/2 /3/2(cA min )e-^ cA — ^ 



(4) 

for the distribution of the smallest eigenvalue. 
Here, J (and I) are (modified) Bessel functions, 
and c = V{nj)ij)). Note that Eq. (||) is simpler 
than the equivalent expression quoted in Ref. Q . 
The RMT result for the connected microscopic 
spectral two-point correlator is 

Ps(x, y) = fd x d y f - d x fd y f , (5) 

where 

/(x,y) = | J dtC(t,2y)-- J dtC(t,2x) 
with 



C(x,y) = 



xJi(x)J {y) - y.J (x)Ji(y) 



x 2 — y 1 



We have performed lattice simulations using — 
2.0 for four different lattice sizes L A with L = 
4, 6, 8, 10. Lattice data for all three quantities are 
compared with the RMT predictions of Eqs. (^)- 
(||) in Fig. |l] for selected lattice sizes. The number 
of configurations is indicated in the figure. The 
agreement is remarkably good even for such mod- 
est lattice sizes. Note that there are no free pa- 
rameters. The only parameter which enters the 
RMT predictions is given by the volume and the 
chiral condensate, which has been fixed by the 
lattice data using the Banks-Casher relation. 

Since the RMT results are obtained in the ther- 
modynamic limit, the quality of the agreement 
between lattice data and RMT increases with in- 
creasing physical volume. Hence, for larger val- 
ues of /3, one requires larger lattices to attain the 
same level of agreement. We have confirmed this 
expectation using (3 = 2.2 on 6 4 and 8 4 lattices 
and (3 = 2.5 on a 16 4 lattice, respectively. 



3. Practical applications 

Having convinced ourselves and, hopefully, the 
reader that the distribution of low-lying eigen- 
values in the limit of large V is given exactly by 
RMT, we now ask how one can make practical use 
of this knowledge. One immediate application is 
the following. We have pointed out earlier that 
the RMT predictions are parameter-free; the only 
parameter appearing in the expressions (essen- 
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Table 1 

The chiral condensate as a function of the lat- 
tice size, determined by the Banks-Casher rela- 
tion (BC) and by a fit to the random-matrix pre- 
diction (RMT). See also Fig. |. 



L 


{H>) (bc) 


(M) (RMT) 


X 2 /dof 


4 


0.1131(19) 


0.1262(17) 


6.75 


G 


0.1209(16) 


0.1263(12) 


1.87 


8 


0.1228(25) 


0.1255(12) 


1.78 


10 


0.1247(22) 


0.1256(10) 


1.15 



tially the chiral condensate) is fixed by the lattice 
data. Given the agreement between lattice data 
and RMT, one can turn the argument around and 
determine the chiral condensate by adjusting the 
energy scale in order to fit the lattice data for, say, 
the microscopic spectral density to the RMT pre- 
diction of Eq. (H). This procedure yields an inde- 
pendent estimate of the chiral condensate which 
should be compared to the value obtained from 
the Banks-Casher relation. This comparison is 
made in Table |] where we also include the % 2 
per degree of freedom for the fit to Eq. (||). The 
same numbers are displayed in Fig. |[ The re- 
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Figure 2. Extrapolation of (4>ip) to the thermo- 
dynamic limit (indicated by the dashed line), via 
the Banks-Casher relation (boxes) and by fitting 
the lattice data for p s (z) to Eq. (|3|) (triangles). 



suits are striking. For all lattice sizes considered, 
the values of (ipip) obtained from the fit to RMT 
agree with each other and with the thermody- 
namic limit within error bars. The values ob- 
tained from the Banks-Casher relation approach 
the thermodynamic limit more slowly. The sur- 
prising result of this analysis is that with a better 
accounting of finite size effects with random ma- 
trix theory, the chiral condensate is remarkably 
volume independent. This approach appears to 
offer an interesting technical advance by making 
it possible to extract certain information about 
the thermodynamic limit from small lattices. One 
can also extract scalar susceptibilities from a sim- 
ilar analysis of spectral two-point functions; work 
in this direction is in progress. Precise numbers 
for the condensate and the susceptibilities as a 
function of temperature are needed to determine 
critical exponents close to the chiral phase transi- 
tion. We are currently testing how the proposed 
approach works at finite temperatures, in partic- 
ular close to the critical temperature. 

4. Summary 

We believe to have demonstrated unambigu- 
ously that the distribution of low-lying eigenval- 
ues of the Dirac operator is described by RMT. 
Some practical applications of this fact have been 
discussed in Sec. |[ We are confident that more 
can be found. 

REFERENCES 

1. H. Leutwyler and A.V. Smilga, Phys. Rev. D 
46, 5607 (1992). 

2. E.V. Shuryak and J.J.M. Verbaarschot, Nucl. 
Phys. A 560, 306 (1993). 

3. T. Banks and A. Casher, Nucl. Phys. B 169, 
103 (1980). 

4. M.E. Berbenni-Bitsch, S. Meyer, A. Sc hafer , 
J.J.M. Verbaarschot, and T. Wettig, hep- 
lat/9704018| . 

5. J.J.M. Verbaarschot, Phys. Rev. Lett. 72, 
2531 (1994). 

6. P.J. Forrester, Nucl. Phys. B 402, 709 (1993). 

7. T. Nagao and P.J. Forrester, Nucl. Phys. B 
435, 401 (1995). 



